Proximity of LaOFeAs to a magnetic instability 
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We investigate the effect of external pressure on the Fe magnetic moment in undoped LaOFeAs 
within the framework of density functional theory and show that this system is close to a magnetic 
instability: The Fe moment is found to drop by nearly a factor of 3 within a pressure range of ± 5 
GPa around the calculated equilibrium volume. While the Fe moments show an unusually strong 
sensitivity to the spin arrangement (type of anti- ferromagnetic structure), the low temperature 
structural distortion is found to have only a minor influence on them. Analysis of the Fermi surface 
topology and nesting features shows that these properties change very little up to pressures of at 
least 10 GPa. We discuss the magnetic instability in terms of the itinerancy of this system. 

PACS numbers: 74.70.-b,74.25.Ha,74.25. Jb 



I. INTRODUCTION 

The discovery of high Tc superconductivity in fluorine 
doped LaOFeAsi with a critical temperature of about 
26 K has stimulated an enormous interest in these com- 
pounds. Shortly after this discovery it became clear that 
a whole family of related compounds shows supercon- 
ductivity at elevated temperatures. Substitution of La 
by other rare earth elements increases Tc up to about 
50 and superconductivity at 38 K was also observed 
in the related Ko.4Bao.6Fe2As2 compound^. 

The undoped parent compound LaOFeAs is an anti- 
ferromagnet with a small ordered Fe moment of about 
0.4 ^b"^- Density functional theory (DFT) calculations 
find, in contrast, a much too large value for the Fe mo- 
ment, close to 2 ^B- The failure of DFT calculations to 
describe the Fe magnetic moment in these systems has 
raised doubts whether the underlying electronic struc- 
ture is correct and whether it provides a sound basis for 
discussion of the superconducting state in the doped com- 
pounds. We argue here that the discrepancies between 
experiment and theory have a physical origin, namely 
the fact that LaOFeAs is close to a magnetic instability. 
Based on our DFT calculations we show that the Fe mo- 
ment is highly susceptible to external pressure and drops 
by almost a factor of 3 within the pressure range from 
-5 to 5 GPa. This drastic change of the Fe moment goes 
along with only subtle changes in the electronic structure, 
what explains the initially apparent differences between 
DFT calculated moments and experimental observations. 
The predicted changes of the Fe moment allow for direct 
experimental verification, either by applying hydrostatic 
pressure or negative pressure which could be realized by 
hydrogenation. 

The crystal structure of LaOFeAs is tetragonal at room 
temperature and consists of FeAs layers separated by 
LaO layers (Fig. [1]). Below Tg ~ 155 K a weak struc- 
tural distortion is observed, followed by the formation 
of a spin density wave (SDW) state around Tjv « 137 
K^. The low temperature crystal structure has been de- 
scribed as either monoclinic^ or orthorhombic^. Both 




FIG. 1: (Color online) Crystal structure of LaOFeAs (a) and 
spin arrangement of the three anti-ferromagnetic structures 
AFl (a), AF2 (b) and AF3 (c) considered in the calculations. 
The unit cell for AF3 is doubled along the c-direction. 



structures differ only marginally from each other, so that 
the symmetry can be described as orthorhombic (see sup- 
plement of Ref. |4j). The anti- ferromagnetic (AF) order 
of most of the undoped parent compounds of the iron 
arsenide superconductors has lead to speculations that 
spin fluctuations could be decisive for the pairing mech- 
anism?. On the other hand, the absence of a SDW state 
has been reported for NdFeAsO, where AF order is only 
observed below 2 K and Fe orders together with the Nd 
moments^. A strong electron-phonon coupling of the Fe 
breathing mode in LaOi_a;Fa;FeAs was reported in Ref.|l, 
which could contribute to the high Tc- Thus, the pair- 
ing mechanism is still under debate and requires further 
experimental and theoretical studies?. 

The effect of hydrostatic pressure on the magnetic 
and superconducting properties of the iron arsenide 
compounds has been extensively studied by experi- 
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mentsiaiiiii^ii^iiiii^ii^iil. For fluorine doped LaOFeAs, 
an increase of the superconducting under pressure 
with a maximum value of 43 K around 4 GPa was re- 
ported-S. For undoped ^Fe2As2 compounds {A—Ga., Sr, 
Ba) , which order anti-ferromagnetically at ambient pres- 
sure, pressure induced superconductivity up to Tc ~ 29 
K for ^=Sr, Bai^ and up to « 12 K for A=Cd^^ 
was observed. Very recently, pressure induced supercon- 
ductivity was also reported for undoped LaOFeAs with 
a maximum value of T„ « 21 K around 12 GPa 
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II. COMPUTATIONAL DETAILS 
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We performed electronic structure calculations in the 
framework of DFT using two high precision all electron 
codes, the full potential local orbital (FPLO) methodic 
and the FLAPW method implemented in WIEN2ki2.. To 
ensure that our conclusions do not depend on the choice 
of functional approximation to DFT, we employed both 
the local spin density approximation (LSDA) and the 
generalized gradient approximation (GGA) . Calculations 
were done for different volumes in the tetragonal crystal 
structure (P4/nmm), as well as in the orthorhombic crys- 
tal structure (Cmma), which are the crystal structures 
above and below T^, respectively. 

We considered three different types of AF spin ar- 
rangements (Fig. [T]): The first cell (AFl) corresponds 
to a checkerboard arrangement in the original unit 
cell, where nearest neighbor Fe atoms are aligned anti- 
ferromagnetically in the xy-plane with a ferromagnetic 
(FM) stacking along the c-axis. Second, we considered a 
stripe-like spin arrangement in the plane with FM stack- 
ing along the c-axis (AF2) in a \/2 x y/2 x 1 supercell. 
The third spin arrangement (AF3) has the same stripe- 
like stacking in the plane as AF2, but in addition the 
spins are also arranged anti-ferromagnetically along the 
c-axis in a \/2 x \/2 x 2 supercell. The experimentally 
observed spin arrangement^ is AF3 and has so far not 
been addressed by electronic structure calculations. 

For all calculations, the scalar relativistic approxima- 
tion was used. The FPLO calculations (FPLO version 
7.00-28) were performed in the local spin density approx- 
imation (LSDA) in the parameterization of Perdew and 
Wang^". For the k-space integrations 512 k-points in 
the full Brillouin zone (FBZ) were used for the struc- 
ture optimization, and the convergence of the magnetic 
moments and Fermi surface properties was checked with 
up to 32768 k-points in the FBZ. In the FLAPW cal- 
culationsi^ the exchange-correlation functional is eval- 
uated within the generalized gradient approximation 
(GGA) , using the Perdew-Burke-Ernzerhof parameteri- 
zation^l. The muffin-tin-radii for La, Fe, As and O were 
chosen as 2.3, 2.15, 2.10 and 1.75 Bohr radii, respec- 
tively. Self-consistent calculations employed a grid of 
4000 (AFl) and 2000 (AF2 and AF3) k-points in the 
FBZ. Rmt X fcmax = 7 was used as plane wave cut-off. 



FIG. 2: (Color online) Left: Fe moment as a function of vol- 
ume (top: LSDA, bottom: GGA) with c/a, zas and ZLa fixed 
to their experimental ambient pressure values for different 
types of AF spin arrangement (see text) . Dashed vertical lines 
denote the calculated LSDA equilibrium volume, the experi- 
mental volume and the calculated GGA equilibrium volume 
(from left to right). Right: Corresponding total energies. 



III. RESULTS 
A. Fe moment as a function of volume 

As a first step, we consider the magnetic properties 
of LaOFeAs as a function of volume in the tetragonal 
structure, with the free parameters c/a, 2;as and ZLa 
fixed initially to their experimental ambient pressure val- 
ues'*. This allows us to distinguish between the effects of 
different spin arrangements, structural parameters and 
the infiuence of the exchange and correlation functional. 
Early electronic structure calculations (see Ref. 22 for an 
overview) showed a confusing variety of results for the 
magnetic properties of LaOFeAs, with values for the cal- 
culated Fe moment between almost zero and 2.6 /is. It 
soon turned out that the Fe moment is highly sensitive 
to the functional as well as the details of the structure 
and the spin arrangement used in the calculations2Si2^. 
While calculations assuming a FM alignment of the spins 
yield magnetic moments of « 0.3 /is (almost in acciden- 
tal coincidence with experiment), calculations using the 
correct AF spin arrangement obtain Fe moments sub- 
stantially larger than the measured «0.4 /i^. 

The variation of the Fe moment and the total energy as 
a function of volume for the three different spin arrange- 
ments is shown in Fig. O The calculated Fe moment for 
the stripe-like spin arrangements AF2 and AF3 at the ex- 
perimental lattice parameters is 1.87 /is within LSDA. 
The corresponding magnetic stabilization energy with re- 
spect to the nonmagnetic state is 3.2 mHartree/Fe, in 
good agreement with the results obtained by Mazin et 
al.^". In agreement with experiment, the stripe-like spin 
arrangement is lowest in energy. Since the magnetic 
coupling between different Fe layers along the c-axis is 
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weak, AF2 and AF3 are very close in energy and have 
also a similar electronic structure and magnetic moments, 
which justifies the use of the AF2 structure in earlier cal- 
culations. 

Although the three AF structures have similar Fe mo- 
ments at the experimental lattice parameters, their be- 
havior with respect to small changes of the volume is re- 
markably different: While the Fe moment of the checker- 
board arrangement AFl sharply drops with decreas- 
ing volume and vanishes already close to the calculated 
LSDA equilibrium volume, it decreases more smoothly 
for the stripe- like arrangements AF2 and AF3. The GGA 
calculations (Fig. [21 bottom) show qualitatively the same 
behavior, which rules out that the observed behavior is 
due to a special property of a certain functional. Within 
GGA, the Fe moments and magnetic stabilization ener- 
gies (2.11 and 6.6 mHartree at the experimental lat- 
tice parameters) are higher than in LSDA, as expected 
from the well-known tendency of GGA to overestimate 
magnetic interactions. 

The remarkable sensitivity of the Fe magnetic moments 
on the type of AF spin arrangement and structural details 
shows already that LaOFeAs is on the verge to a magnetic 
instability with respect to changes of the volume. In the 
vicinity of a magnetic transition or instability, parameter- 
free LSDA or GGA calculations can not be expected to 
yield the exact value for the magnetic moment, but they 
should qualitatively reproduce the behavior as a function 
of an external parameter like pressure. However, the re- 
sults shown so far do not yet fully explain the deviations 
between the calculated and measured Fe moment: While 
a sharp drop as in the case of AFl could well explain 
those deviations, the Fe moments remain substantially 
too large in the vicinity of the equilibrium volume for 
the correct AF3 spin structure. 



B. Structural optimization: Fe moment as a 
function of pressure 

Our next step is to consider structural optimizations 
of the free parameters c/a, zas and ZLa for different vol- 
umes. In earlier publications it was pointed out that the 
magnetic moments in LaOFeAs are highly sensitive to 
structural parameters, especially to the height of the As 
position ZAsT^^- Under pressure, changes of the struc- 
tural parameters with respect to their ambient pressure 
values can be expected. In the following we restrict our- 
selves to LSDA calculations, which are better suited to 
describe the magnetic behavior in the vicinity of a mag- 
netic instability, due to the tendency of GGA to overesti- 
mate the magnetic interactions. LSDA calculations have 
for example successfully been used to predict a meta- 
magnetic transition in YC05 under pressure^^. 

The structural optimizations have been performed 
with spin-polarized calculations in the AF3 structure. 
The c/a-ratio shrinks considerably with pressure, reflect- 
ing a weak inter-layer coupling, and also the internal 
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FIG. 3: Fe moment as a function of pressure for the AF3 spin- 
structure with optimized structural parameters. The inset 
shows the variation of c/a, V, zab and Zha. with pressure. 

parameters are subject to considerable changes (inset 
of Fig. [3]). The Fe moment, calculated with the opti- 
mized parameters as function of pressure, is shown in 
Fig [3l Within ± 5 GPa around the calculated equilib- 
rium volume (zero pressure), the Fe moment drops by 
nearly a factor of 3. The calculated moment at zero 
pressure is 0.75 /is and thus still about two times larger 
than the experimental value. At a pressure of about 5 
GPa, the calculated Fe moment coincides with the one 
observed in experiments. The magnetic stabilization en- 
ergy decreases from 2.6 mHartree/Fe at -10 GPa to about 
0.3 mHartree/Fe at ambient pressure. Spin fluctuations, 
which are only incompletely included in LSDA or GGA 
calculations, are expected to lead to a substantial sup- 
pression of the magnetic moment when the magnetic 
stabilization energy is of the order of 0.5 mHartree per 
atom^^. Hence, an even sharper reduction of the Fe mo- 
ment with pressure than the one shown in Fig. [3] might 
be observed in experiment, although the precise effect of 
the spin fluctuations can not be estimated. 

Up to now, we did not consider the effect of the or- 
thorhombic lattice distortion observed at low tempera- 
tures. We find indeed a minimum in the total energy 
with a small deviation in the 6/a-ratio of about 1% from 
a tetragonal lattice, in agreement with experiment and 
an earlier report by Yildirim^i. However, in contrast to 
Yildirim's work we find only a minor influence of this 
orthorhombic distortion on the Fe moments (~ 0.05/1^), 
leaving the data shown in Fig. [3] basically unchanged. 



IV. DISCUSSION 

The result of our study shown in Figs. [2] and [3] pre- 
dicts that LaOFeAs is close to a magnetic instability. At 
a pressure of about 5 GPa, the Fe moment calculated 
within LSDA coincides with the one observed in experi- 



ments. This means that LaOFeAs is on the right side of 
the transition shown in Fig. [31 A strong increase of the Fe 
moment would be expected for negative pressure condi- 
tions, which are of course not straightforward to reahze 
in experiments. However, it is weh known that hydro- 
genation can lead to a sizeable increase of the volume 
and could thus serve as a medium to simulate negative 
pressure for LaOFeAs. Since the electronic structure of 
different iron arsenide compounds is quite similar-^-, we 
would also expect that the behavior we found in LaOFeAs 
can be observed in other iron arsenide compounds. 

While finalizing this work, further studie o^^'^^i'^^ on the 
magnetic properties of iron arsenide compounds under 
pressure were reported, which support our interpreta- 
tion. Kumar et al^ have performed a combined the- 
oretical and experimental investigation for SrFe2 As2 . In 
their LSDA calculations, they find a suppression of the 
magnetism at a critical pressure of about 10 GPa, which 
is slightly higher than the critical pressure extrapolated 
from experiments of 4-5 GPa. Yildirim^ has performed 
GGA pseudopotential calculations for CaFe2As2 and 
LaOFeAs under pressure. His calculations for CaFe2As2 
show a suppression of the Fe magnetic moment around 
10 GPa (note however that CaFe2As2 is special due to 
the presence of a collapsed tetragonal phase). In the case 
of LaOFeAs, the Fe moment calculated within the GGA 
pseudopotential approach is higher at ambient pressure 
than in our LSDA calculations and remains close to 2 /is 
up to 10 GPa, but drops to zero at 20 GPa. However, 
the data shown do not allow to judge if the transition is 
smooth as in our calculations or if a sudden collapse of 
the moments under pressure occurs. Finally, Xie et al.^^ 
have performed calculations for BaFe2As2 under pres- 
sure and found a suppression of the Fe moment around 
13 GPa. 

An explanation of the magnetism in LaOFeAs in terms 
of localized magnetic moments is difficult, if not impossi- 
ble. First, the magnetic moments are soft and depend on 
the spin arrangement and structural details, which is not 
compatible with a simple Heisenberg model. Second, the 
total band width of the Fe 3d states amounts to about 7 
eV. Near the Fermi energy, all five d-orbitals contribute 
to the density of states (DOS), with little admixture of As 
4p states. In localized systems, crystal field splittings are 
a valuable tool to predict the spin state^i^. The related 
crystal field splittings of the Fe 3d states in LaOFeAs 
(evaluated from the center of gravity of the correspond- 
ing partial DOS) are well below 0.5 eV and thus much 
smaller than the band width. 

In an itinerant magnet like LaOFeAs'^'*, the magnetic 
state is determined by a delicate balance between ki- 
netic energy (favoring a nonmagnetic state) and the gain 
in exchange energy by spin polarization. With increas- 
ing pressure, the bands of LaOFeAs are broadened and 
weight is shifted away from the Fermi energy. A (rough) 
quantitative measure for this shift can be obtained from 
the integrated partial DOS weighted with a Gaussian 
around the Fermi energy^^, which yields a 20-30% re- 
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FIG. 4: (Color online) Cuts through the FS of LaOFeAs for 
two different volumes (top: V=141.9 A^, bottom: V=120 
A^). Cuts perpendicular to the c-axis in the F-plane (left) 
and in the ^-plane (right) are shown. To visualize nesting, 
the same cuts shifted by Q = (tt, tt, 0) are also drawn (light 
color). Similar nesting is observed throughout V-Z. 



duction from -10 to +10 GPa, where the Fe moment 
essentially drops from 2/iB to almost zero. We do not 
observe abrupt changes in the electronic structure in this 
pressure range, consistent with the smooth (but rapid) 
decrease of the Fe moment in Fig. [31 However, close to 
the magnetic instability, the observed changes are suffi- 
cient to alter the magnetic state. 

Mazin et al^ have pointed out that the stripe-like spin 
arrangement is stabilized by nesting features in the para- 
magnetic Fermi surface (FS). This is confirmed by our 
calculations, which also show that the related nesting 
features remain robust under pressure. The FS consists 
of five sheets, with two cylindrical hole sheets around F 
and two cylindrical electron sheets around M, nested by 
a vector Q = (tt, tt, 0). In addition, there is a hole pocket 
around Z, whose shape depends strongly on structural 
details'^^ . The larger one of the F-centered FS sheets be- 
comes more three dimensional around 10 GPa, but the 
topology of the Fermi surface does not change up to pres- 
sures of at least 10 GPa, and also the nesting features are 
surprisingly robust (Fig. [4]). The nesting is never perfect, 
but remains substantial throughout the considered pres- 
sure range, which explains the relative stability of the 
magnetism in the AF3 structure. 



V. CONCLUSIONS 

In summary we have shown that LaOFeAs is close to a 
magnetic instability, which explains the discrepancies be- 
tween the values for the Fe moment found in experiment 
and DFT calculations. On the basis of our calculations 
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we expect a strong increase of the Fe moment with in- financial support through the TRR/SFB 49 program, 
creasing volume, which could be realized for example by 
hydrogenation. The Fermi surface topology and the re- 
ported nesting properties are fairly robust up to pressures 
of at least 10 GPa. 
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